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ABSTRACT 

Coronal emission line intensities are commonly used to measure electron tem- 
peratures using emission measure and/or line ratio methods. In the presence of 
systematic errors in atomic excitation calculations and data noise, the informa- 
tion on underlying temperature distributions is fundamentally limited. Increas- 
ing the number of emission lines used does not necessarily improve the ability to 
discriminate between different kinds of temperature distributions. 

Subject headings: Sun: corona 



Introduction 



Spectroscop ic measurements of the t emperature of c orona l plasma have been made 
for decades fe.g. ISeaton Il962l : iNocilhoO.sk iFeldman et all Jl999h : iFeldman fc Landil fhoosl ) 
and Landi fc Feldman ( 20081 ) presented intriguing evidence that electron temperatures of 
solar coronal plasma, measured from emission line spectra obtained high in the coroma, 
are clumped into several peaks, and are not broadly distributed. Narrow distributions of 
plasma temperature have important implications for the energy balance of the corona. In 
the methods used, frequency integrated line intensities / are assumed to be simple func- 
tions G(T) of the logarithm of the electron temperature, T, because of the dominance of 
two body collisional processes whic h lead to the well-known "coronal approximation" (e.g. 
Woolley fc Allen! Il948l ; ISeatonlll964l ). Two approaches were used. In one they solved for the 
differential emission measure which is an optimal solution to the inverse problem: 



Gi(T)£(T)dT, i = l...n 



1 The National Center for Atmospheric Research is sponsored by the National Science Foundation 



-2- 



where there is a set of n different emission lines. The other method sought the single 
logarithmic temperature T such that £(T) oc 5{T — T ), by plotting Ii/Gi(T) and identifying 
i ntersection points for the observed lines. These approaches are, in fact, formally equivalent 



( IMcIntosh et al.lll998l ). In both cases, information must be added to obtai n both £(T) and 



T n : in the first case the problem has to be "regularized" as it is ill-posed (e.g. I Craig fc Brown 



19861 ). one searches for example for the "least structured" solution £(T) that is compatible 



with the data. In the second case it is assumed that t he plasma is indeed approximately 



isothermal. In the analysis of iFeldman fc Landil (120081 ). the individual peaks in the £(T) 



functions have widths of 0.1-0.2 in T, and the curves of Ii/Gi(T) intersect one another 
within similar margins. 

Some general questions arise. Given typical uncertainties in observable and model pa- 
rameters, are the observed data compatible with different £(T) distributions? What, then, 
is the accuracy of the derived temperatures? In this paper these questions are addressed 
by asking, how broad can £(T) functions be to be incompatible with observed data. By 
how much can one change the temperature of an isothermal plasma before the differences in 
line intensities become significant? In the problem at hand there are unavoidably large and 
systematic uncertainties in the G(T) functions. These uncertainties limit the information 
which can be extracted from emission line spectra. It is found that acceptable widths w of 
the £(T) functions exceed the precision by which the lines can in principle determine that 
two different plasmas have slightly different temperatures. Thus, these widths set the lower 
limit to the ability of emission line techniques to diagnose electron temperatures. 



2. Calculations of temperature sensitive lines 
2.1. Uncertainties in G(T) 



G(T) functions depend linearly on electron impact excitation rates, on elemental abun- 
dances, and on factors influencing the ionization fraction for a given line. For coronal ions 
uncertainties in electron impact excitation cross sections, based primarily u pon sophisti- 



cated calculations, are typically ±10% but can be considerably higher (e.g. IStorey et al. 
20021 ). Other uncertainties in G(T) arise from element abundance uncertainties, bound- 
free cross sections (ionization, recombination), and can include systematic errors from non- 
ionization equilibrium effects, from any significant non-thermal populations of electrons, 
modifications of dielectronic recombina tion rates due to finite density plasmas (e.g. lSummers 
1974] ; iBadnell et al.l 12003c I Judge! 120071 ). and even radiative transfer. Given these consider- 
ations, we adopt fractional uncertainties of e = 0.2 for all G(T) functions, which should 
probably be considered a lower limit. 
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2.2. Numerical calculations 

C alculat ions were performed for lines of Fe VIII - Fe XVII using the DIPER package 



( |Judgdl20071 ). The upper left panel of Figure [T] shows G(T) functions computed for a typical 
coronal electron pressure of 0.1 dyne cm -2 . In constructing this figure, intensities Ii SO th were 
computed using an isothermal plasma at 1.3 x 10 6 K. Intensities were then recomputed using 
a Gaussian function with unit (— oo — > +oo) integral: 

m = e-(^) (2) 

y'TTW 

The intensities were recomputed, and a value of \ 2 computed using the above uncertainties 
in G(T) added in quadrature to the uncertainty from photon counting statistics, assuming 
that the brightest line in the spectrum has accumulated 10 4 counts. These uncertainties 
are not unrealistic values for what might be expected from UV and EUV spectrographs of 
the modern era. This particular calculation shows that an isothermal plasma with T = T 
cannot formally be distinguished from a plasma with Gaussian distributions of width w less 
than 0.1. Informally, the figure shows a systematic increase in the ratio I /hsoth in lines with 
G(T) peaking away from To. One may be tempted to claim that the curvature indicates 
the presence of broader distributions. But these calculations are made with the same set 
of G(T) functions, where in reality one would compare these theoretical values with those 
(unknown) functions, free of systematic errors, present in the real Sun. 



2.3. Analytical approximations 

To generalize these results, analytical results can be obtained using suitably simple 
functions to approximate Gj(T). Here, Gaussian functions are adopted both for £(T), as in 
eq. ©, and G^T): 

( T-Ts \ 2 

Gi(T) = a t e~y—) . (3) 

Very similar results are obtained using £(T) represented by a "top-hat" function of full 
width 2w. These were chosen because they are not qualitatively dissimilar from the shapes 
of typical lines, and they are analytically simple. Then 

OO /t_t. \2 - T _ T()X 2 



Ii = -7^/ e \ «* ) e-y-^) dT, i = l...n (4) 

y TTW ' 



oo 



y^T 
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When w — > 0, then Jj = If = a^e \ > , i.e. the isothermal case. If w 2 <C v\ and for 
source temperatures Ti = T + r (r <C |T — T|- large values of |T — T| are needed for 
significant temperature sensitivity between at least some lines), these expressions can be 
expanded as a Taylor series, to yield changes in Ij as functions of small values of w and S. 
The first derivatives are 

I = 2 ^ <«> 
^ = 27° W 1 P»- T °) 2 ^ (7 ) 

Dividing by Ii, the first order fractional changes in 7j resulting from small changes r and w 



are: 



1 dh „ T - T 



1 S/j 2u> 2 f (Ti — Tp) 2 1\ 

A - = 7F^ "T (— 2) (9 > 

These correction terms can be used to compute the reduced x 2 statistic, to first order, using 
the uncertainties e,-: 



X n - 1 ^ Ihl n - 1 ^ e? 



-i 



Define a;, = |T — Ti\/vi, which measures the difference in temperature of the peak of the 
source relative to the peak of line i in units of the Gaussian width of G(T). Any sensible 
set of lines must include values of Xi varying from to numbers ^> 1. But for xi > 3, the 
lines become very weak with accompanying large observational uncertainties. Practically 
speaking, we can consider most values of Xi to be spread between and 3. Any given set of 
lines cannot be used to discriminate an isothermal source at temperature To from another 
at temperature T ± r when (using A = A T ), the majority of lines satisfy 

\A < 7^- Vi, V % = 1 . . .71. (11) 

2xi 

Similarly one cannot discriminate an isothermal source at temperature To from a broadened 
distribution with Gaussian parameter r when (using A = A w ) 



w ~ v \ 2x2 ~ x l Vu v i = 1 ' ' ' n ' ^ 

Using Xi ~ 1, Vi ~ 0.15, e ~ 0.2, the set of emission lines lines has sensitivities to the peak 
of the temperature and its width given by w of 

r~ 0.015, w~0.13. (13) 
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This result for w is in reasonable agreement with the numerical calculation. The "bluntness 
of the thermometer" is, in this approximation, measured by w. This quantity therefore 
determines the accuracy of temperature measurements. The accuracy is a factor of 8 worse 
than the precision with which different isothermal plasma temperatures can be discerned, 
relative to one another, which is represented by r. The width w is weakly dependent on 
the intrinsic uncertainties £j, linearly dependent on widths Vi of the Gi(T) function, and 
inversely proportional to 

Given the weakness of lines with large Xi values, the nature of the uncertainties it 
relatively little can be done to increase the thermometer's accuracy. However, one can try 
to include very strong lines (large a^) which are formed far from their peak temperature 7$, 
i.e. when x-i > 3, say, such that the counts are still sufficiently large to make q approach 
the systematic uncertainties in atomic calculations. One example might be the strong line 
of Fe IX near 17.1 nm whose G(T) function is prominent in Figured! In such cases 



Note however, that such lines tend also to have large widths v^. 

Temperature estimates cannot however be improved simply by analyzing more lines 
with similar values of Xi, because then the largest acceptable width w is almost independent 
of the number of lines n used. The uncertainties in q are, in an important sense, irreducible, 
as, even in the limit of perfect observations, they result from uncertainties in atomic cross 
sections and in our lack of knowledge of the conditions controlling level populations in the 
solar coronal plasmas. 



The largest values of w compatible with observations yield the accuracy with which 
emission lines can measure electron temperatures. The widths of the distributions compatible 
with known sources of uncertainties set a natural limit on the sharpness of a detectable peak 
in the underlying emission measure distributions. Typically the full width, 2w, is 0.2 to 
0.3 in the logarithmic electron temperature. The precision is an order of magnitude better, 
being estimated using the sensitivity of the spectra to two isothermal plasmas differing in 
logarithmic temperature by r. Although this implies that two strictly isothermal plasmas 
could be differentiated by using line ratios to a precision of order 6^/2 ~ 0.015 in the 
logarithmic electron temperature, we will perhaps never know if such isothermal plasmas 
exist with a width narrower than to, based upon these data. 





(14) 



3. 



Discussion 



These properties prompt the following comments: 
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Error bars quoted below ~ 0.1 in logarithmic temperatures are not credible using these 
techniques. Landi and Feldman quote error bars of 0.04 and 0.05 in their work which 
seeks to provide spectroscopic evidence that the coronal temperature is "quantized" . 



2. The widths (FWHM) of peaks in the f (T) functions found bv lLandi fc Feldmanl (12008 
their figure 5) are between 0.1 and 0.2. These widths are characteristic of the limit 
with which emission lines can determine the isothermality of the emitting plasma, 
and are therefore determined more by regularization than by the data themselves. 
Their analysis is consistent with isothermal plasma, but is limited by the uncertainties 
discussed here. 

3. Coronal temperatures are controlled by (unknown) h eating mechan i sms a nd cooling 



by heat conduction, flows, radiation. The scalings of iRosner et al.l (119781 ). based on 
simple energy balance considerations, show that a given coronal electron temperature 
T e in a loop of length L requires an energy flux density 



T °c (T e L) 



7/2 



to sustain it agains conductive and radiation losses. Therefore, for a given L, if T e is 
uncertain to ±0.15 in its logarithm, the absolute value of logio^ 7 can be determined 
from a given length loop to an accuracy of ~ ±0.5. The energy flux is constrained only 
to somewhere within a range spanning a factor of ten using emission lines. Conversely, 
variations of the energy flux T dissipated in coronal loops by < 0.5 in the logarithm 
cannot be distinguished through temperature dependent emission line methods. 

4. Relative changes in T could, however, be detected to within about ±12%, but again 
this precision depends on the plasma being nearly isothermal (i.e. actual w ~ t), 
which is neither known from the data nor expected from first principles. 



The analytical results offer some hope that carefully selected sets of emission lines can 
be used to reduce the widths of distributions compatible with the data. Specifically, using 
lines within the same element removes errors arising from incorrect abundances, and can 
include lines with large values of \xi\. In this case just a few lines with different values 
of Xi can reduce w via equation (fT4j) ). Using lines from within the same ion stage yields 
significantly smaller values of ~ 0.1, but while the precision improves linearly with e» the 
accuracy improves only quadratically. 



The author is grateful to Scott Mcintosh for carefully reading the manuscript. 
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Fig. 1. — G(T) functions for Fe VIII- Fe XVII are shown as a function of T (logarithm of 
the electron temperature) in the upper left panel. Solid and dashed lines mark ions with 
odd and even charges repectively, for clarity. The other panels plot the ratios of intensities 
computed using the Gaussian functions shown (solid lines) divided by the intensities for an 
isothermal calculation centered at the same temperature (in this case T =log 10 1.2 x 10 6 
K). The uncertainties are set to > 20% of the computed line intensities (see text) and the 
computed x 2 parameters are listed with the Gaussian width w. 



